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Abstract 

We study a system of N Bose atoms trapped by a symmetric harmonic 
potential, interacting via weak central forces. Considering the ground state 
of the rotating system as a function of the two conserved quantities, the total 
angular momentum and its collective component, we develop an algebraic 
approach to derive exact wave functions and energies of these ground states. 
We describe a broad class of the interactions for which these results are valid. 
This universality class is defined by simple integral condition on the potential. 
Most of the potentials of practical interest which have pronounced repulsive 
component belong to this universality class. 
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I. INTRODUCTION 



Recent progress in creation of magnetic traps for cold atoms made possible the study of 
the effects of quantum degeneracy in systems of finite numbers of interacting identical par- 
ticles 0. One of the most interesting phenomena is the Bose-Einstein condensation, which 
can be studied in different coupling regimes in such dilute atomic gas systems. Theoretical 
studies of such systems of weakly interacting bosons, confined by parabolic potential, have 
become very active 

Response of such system to rotation and the onset of vorticity in the condensate are 
among the most interesting questions attracting attention of both experimentalists and 
theorists @-^|. In this connection, it is very important to study the structure and the 
spectrum of the ground states of rotating condensates at a given angular momentum, called 
the yrast states see Fig. 1. The notion of yrast states is borrowed from nuclear physics, 
where the states of many-body system having highest spins at a given excitation energy play 
also a special role [TIJ . 



A particular problem of wide recent interest arises in the so-called weak coupling 
limit, when the repulsive interactions between bosonic atoms can be considered weak as 
compared to the spacing between different oscillator levels in the harmonic trap. This limit 
is expected to be reached in future experiments. In particular, an interesting possibility 
to study different coupling regimes is related to using the Feshbach resonance |lT| , [l2l to vary 
the strength of the effective interatomic interactions. From the theoretical viewpoint, the 
advantage of the weak coupling limit is that it allows analytical treatment of the problem. 

We will consider here N spinless Bose atoms in a spherically symmetric harmonic trap, 
assuming the interactions to be weak. It is sufficient to consider the two-dimensional case, 
because the three-dimensional case can be reduced to the former in the weak coupling limit, 
as was shown in Refs. @||. It is expedient to start with the case of noninteracting particles. 
The noninteracting system has equidistant spectrum Tium where u is the trapping oscillator 
frequency and n an integer. Each level is degenerate, and the degeneracy grows exponentially 
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with n for n — N ^> 1 jp. This degeneracy is related to the number of ways to distribute 
the total energy hurt among the Bose particles. The short-range interactions V(r) between 
the atoms are assumed weak in the sense that hoppings between different hurt levels can be 
neglected, 

N(V) < hu, (1) 

where (V) is a typical matrix element of the interaction [|J]. The problem is therefore to find 
a nonperturbative solution for the highly degenerate states at a single level hurt, which is 
similar in spirit to the problem of the lowest Landau level for the electrons in high magnetic 



field, which arises in the theory of the fractional quantum Hall effect fllqjbj or to the problem 
of compound states in an atomic nucleus. The yrast states are those with minimal energy 
at given angular momentum, L. This is illustrated in Fig. 2, where the spectrum patterns 
are shown for the cases of zero and finite interaction. 

As is usually the case for interacting many-body systems, the evaluation of the exact 
ground state is a prohibitive task, even with the simplification introduced by the weak 
coupling limit []Hj-^5[. The yrast states in the case of attractive <5-forces have been 

found analytically by Wilkin et al. ||. Later, these results have been shown valid for 
a broad class of attractive interactions in Ref. [p5 l. The case of repulsive interaction is 



more difficult to analyze P-[^|i5H2l|j . One of the first important results obtained for the 



repulsive forces was that of Bertsch and Papenbrock ||: these authors diagonalized the 
repulsive (^-interaction numerically and suggested analytical formulae for the wave functions 
and energies of the yrast states. Later, it was shown analytically by several authors using 
various methods [18[ , ^6| , |17 , ^5) , 19[ , that the states of the form of Bertsch and Papenbrock are 



indeed eigenstates of the Hamiltonian. Only recently, it was shown that these states indeed 
correspond to minimum energy P2| , ||23|| . 

In this work, we provide rigorous analytical solution for the yrast states. In fact, we 
consider a more general problem []22|, of finding the ground state as a function of two 



quantum numbers, the total angular momentum, L, and the angular momentum of internal 



excitations (we discuss this quantity in the following sections). Such "generalized yrast 
states" include the usual yrast states as a limiting case. These solutions are valid in the 
region L < N . Physical interpretation of the results is quite transparent. In fact, the form 
of the yrast wave functions which was drawn from numerics || for the case of ^-interaction, 
turns out to be valid for a broad class of repulsive interactions. The universality class of 
such interactions [which includes, but is not limited to, delta-forces, Gaussian forces with 
arbitrary range, Coulomb (1/r) and log-Coulomb [log(r)\ forces] is described by an explicit 
sufficiency condition. 

The structure of the paper is the following. In section II, we present detailed formulation 
of the problem in the weak coupling limit, introducing the relevant "partition subspace" and 
the basic ingredients of the following consideration. In section III, we derive the operator 
expansion for the Hamiltonian projected onto the partition space. In section IV, we dis- 
cuss the additional conserved quantum number, "seniority", which can be interpreted as a 
collective contribution to the angular momentum, and introduce the notion of generalized 
yrast states. The method of algebraic decomposition of the Hamiltonian is introduced and 
discussed in section V. Its application to the present problem is developed in section VI. 
Section VII is devoted to study the properties of eigenvalues of the "perturbation" defined 
in the section VI. The results for the generalized yrast states are presented, in general form, 
in section VIII. In section IV, we consider various examples of the interaction potentials and 
discuss the results. Section X is devoted to detailed discussion of the applicability condition 
and of the corresponding universality class of attractive potentials. Section XI summarizes 
the work. 

II. THE PARTITION SUBSPACE 

In the weak coupling limit discussed in the previous section, the problem splits onto 
series of independent problems for each value of the total angular momentum L. At given 
angular momentum L, the Hamiltonian (H = m = 1) is the sum 
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H = luH + V, (2) 

where the first term (the c-number equal to the energy of degenerate level) comes from the 
noninteracting Hamiltonian in the parabolic trap, 

N / -Q 2 \ 

^o = E(| + Y^J - ^o, (3) 

with pi and r*j denoting the momentum and coordinate of the i-th particle. Hereafter, we 
set oj = 1 for brevity. The second, nontrivial, term V stems from the two-body interaction 

V = J2V( rij ) -> V, (4) 

i>j 

projected onto the single-level problem, which is worked out below. 

In two dimensional problems involving harmonic potentials, it is convenient to use the 
complex variable z — x + iy and the conjugated z* = x — iy instead of the position vector 
r — (x,y). Introducing the notations = §(J|— i-§^) and = |(Jj + ? J^)> it is convenient 
to employ the tetrad of ladder operators a + , a, b + and b 

, z d ,j. 2* 9 



2 9^*' 2 
a=(a + )t = f + |, 6^ + )t.| + JL (5, 

for each particle. For brevity, the particle markers are suppressed in Eqs.(|5|), dagger denotes 
Hermitean conjugation. The role of these operator is to raise (lower) the powers of z% and 
z* in the preexponentials of many-body wave functions, which are all polynomials times the 
Gaussian factor |0) = exp(— 1/2 \ z k\ 2 )- Thus, we have 

Zl \0) = at\0), $\0) = b+\0), 

a,i\Q) = 0, bi\0) = 0, etc. (6) 
The only nonzero commutators between (|5|) are given by 

[ai,a+] = [bi,bj] = (7) 
In two dimensions, the total angular momentum is the difference 



L = L + - L_ (8) 

of numbers of "up" and "down" quanta, L + = J2i a t a k an d L_ = J2i H^k, respectively. 
At the same time, the noninteracting Hamiltonian on the left hand side of (|3p is given by 

H = L + + L„+ N. (9) 

By definition, the yrast states have minimum energy at given L. They must therefore belong 
to the subspace with 

L_ = 0, L + = L. (10) 

The first term in (|||) is therefore reduced to a constant given by 

H = N + L. (11) 

The subspace defined by (|T0l) is spanned by the homogeneous symmetric polynomials of 
degree L, polyg(zi) which do not involve conjugated z*'s at all, 

m{L)=poly L s {atm (12) 



The orthogonal basis of these symmetric polynomials can be constructed as follows |26 
Consider the monomials 

m = {h,l 2 ,...,l N } = 4z l i...z 1 ^ (13) 

corresponding to partition of integer L, J2^n = L. The states m|0) with different sets 

n 

{li, I2, In} are mutually orthogonal. The basis symmetric polynomials are given by 
symmetrized linear combinations of flT3| ) 

[h,l 2 ,...,l N ]=P s m, (14) 

denoted by corresponding ordered partition 1 2 , In] with l x > l 2 > ... > In- Here, the 
operator of symmetrization P$ is the standard sum over all the permutations. In the case 
L < N which we are interested in here, the number of the basis states ([14]) is given by "the 
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number of unrestricted partitions" [h,h, ■■■,In] of positive integer L p7j . At L < N, this 
number does not depend on N; hereafter, it is denoted by p(L). For example, at iV = L = 3 
the complete basis ( |H| ) of the partition space (|T2| ) is spanned by the following states 

[1,1,1] =Z!Z2Z 3 = 10,3,0,0...}, 

[2, 1, 0] = z\z 2 + z x z\ + z%z 3 + z 2 z\ + z\z\ + z 3 z\ = |1, 1, 1, 0...), 
[3, 0, 0] = z\ + z\ + zl = |2,0,0,1...), (15) 

which correspond to the number of partitions p(3) = 3 (normalization factors and the 
Gaussian, |0), are suppressed for brevity). In Eqs.(|l|), the right hand side of each line 
gives interpretation of many-body state in terms of the single-particle occupation numbers 
of the oscillator states |no, ni, ...) with rtj the number of Bose particles in the state with 
i oscillator quanta. 

The numbers of partitions, p(L), can be computed from the generating function, using 
the formula 

oo oo -i 

L=0 k=l L 1 



An explicit (however involved) expression for p(L) is given in p7|. For small L, the values 
of p(L) are 

p(2) = 2, p(3) = 3, p(4) = 5, p(5) = 7, 

p(6) = ll, p(7) = 15, p(8) = 22, ... (16) 



At L large, the number of states grows exponentially |27 



p(L>l)oce L ^, (17) 



as is very typical for many-body systems JT0| 



The problem of finding the spectrum consists of diagonalization of the interaction, the 
second term in (0), in the "partition subspace" ([T2]Jl^) ; the eigenvalues of (@) have the form 

S = N + L + E } (18) 



where E is the interaction contribution. It should be noted that this term can not be found 
with any kind of perturbation theory, as it comes from diagonalization of the interaction 
within the subspace of degenerate states. 

In this paper, we restrict ourselves to the case L < N only. At L > N, the whole situation 
is changed drastically. First, the number of basis states is not given by p(L) anymore and it 
depends on both L and N. For example, the simplest states with partition [1, 1, 1...1], which 
play the role of generating functions for the ground states at L < N (see section V), can not 
be constructed at L > N. It is therefore reasonable to expect something similar to "phase 
transition" at L = N. In the earlier studies ||, the signatures of this phase transition have 
been observed in numerical simulations. 



III. THE INTERACTION WITHIN THE PARTITION SUBSPACE. "OPERATOR 

EXPANSION" 

In this section, we work out a convenient representation for the V, Cf. (^), the interaction 
V projected onto the "partition space" which was described in the previous section. Using 
condition Eq. fllOl) , the projected interaction V in Eq.(|2|) can be written |T5lj22 



as 



V= E PoP£VP+P ~, (19) 

L=— oo 

where we use the standard number-of-quanta projector 

P£ = — dae ia ( L -^< a >>) , (20) 
2tt jo 

and the similar expression for P£ with substitutions a~l — > bf and — > hi. Due to construc- 
tion, the transitions between different L-sectors are suppressed in (0), and the resulting 
effective interaction (^) applies equally to every L-sector. 
With the help of (]2TiD, Eq.([TP|) can be cast in the form 

oo 

V = — I daP^e ia ^i a t ai Ve- ia ^i a t ai Po. (21) 
2n J 



-'TO 



s 



In order to evaluate (pi]), we use the Fourier representation of the interaction in the right 
hand side of (f21|), 



oo oo 



oo — oo 

where 



V = Y, J <k* J dqye^-^+^+o+^+^Wq (22) 



Vq = 2~?r Jo TdrJ ° ( r ^J q * + q y) V ^ 



is the two-dimensional Fourier transform of the central potential V(r). Here, Jo is the 
Bessel function |27j and notation q± = (q x ±iq y ) / \/2 is introduced. The two-particle operator 



combinations 

1 , 1 1 _ i 1 



4 = bt^-^ibt-bi), 



a± = (a±y, b± = (b±y (23) 



came from resolving x,,^ and Xj,yj from (g). Substituting fl22|) in flT9|), we make use of 
Baker-Hausdorf relation 



e a+-a = e L_J ga+e _ a 



which is valid for any pair of boson operators a + and a. Expanding ( ^2|) to powers of q 2 and 
evaluating integrals term by term, we use the relation 

/•oo 

/ dqe" 3 /2 q 2n+1 J (qr) = 2 n n\M(n + 1, 1, -r 2 /2) 
Jo 



where M(/i, i/, x) is the Kummer confluent hypergeometric function |27]. Proceeding in this 
manner, we obtain expansion of V 

2[L/2] 



V 



E (~l) k s k B k = N(N 1} - SxS 1 + ... + (24) 



in terms of the two-particle normal-ordered operators defined as 

JV 



£ fc = E4> ''U "U"'ir ( 25 ) 

j>i 



In Eg. (|24|) , [q] denotes integer part of real number q. The "strength parameters", s&, are 
related to the potential, V(r), via the integrals 



1 r°° 

s k = -n dtM(k + 1,1, -t)V{V2t). 
k\ Jo 



(26) 

with M the Kummer function ||27|| . For k integer, the Kummer function can be expressed 
in terms of the Laguerre polynomials L° 

M(k + 1,1, -t) =e~ t Ll(t). 



In the expansion fl24]) , the highest possible order k of the operators B k is L for L even, 
and L — 1 for L odd. In the lowest order term, no operators are involved, and we have 
B° = JV(JV-l)/2. 

Eq.(|24]) can be written in another convenient form 



[L/2] 

v=J2v h , 

k=0 

V k = s 2k (B 2k - B 2 '- 1 ) + (s 2k+2 - s 2k+1 )B 2k+1 . 



(27) 



Here, = 0. This second form of the expansion contains operator structures collected 
in the way convenient to apply the method of "algebraic decomposition" described in the 
following sections. 

The operator expansion (24j27j) is universal, while the particular shape of the potential 



V(r) is described by the integrals with Kummer function M [27]. We should stress that 



expansion (E7h is exact for any interaction V(r) whose moments s k are finite [15,24 . 



IV. THE QUANTUM NUMBER "SENIORITY" , GENERALIZED YRAST 
STATES AND CORRESPONDENCE RULES 

In this section, we discuss an extra conserved quantum number the system enjoys, in 
addition to the energy and the total angular momentum. @,f|]. This quantity, which is very 
important for classification of states, is reminiscent to the "center-of mass" mode in nuclear 
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problems [10]. Here, it can be interpreted as collective contribution to the total angular 



momentum L. Indeed, the pair of the mutually conjugated collective ladder operators 

N a + N a- 

A+ = J2^ A = y-^ [ AA +] = 1 (28) 



commute with any two-body combinations in (]27|) as well as with the angular momentum, 

[A, a±] = 0, [A,L + ] = 0, [A,L_] = 0, (29) 

and they therefore commute with both terms in the Hamiltonian H (|2|), thus we have 
[A, H] = 0. The number of collective quanta, A + A, is therefore a conserved quantity. The 
mutual eigenfunctions of the triad of operators H, L and A + A = v can be found in following 
factorized form 

**(L,v)= Z poly L s ~ v {~ Zl ) e-E^, (30) 
where Zi = Zi — Z and 



A? 



Z = J2*i/N (31) 

i=l 

is the collective variable. The additional index k stands to distinguish between the different 
states in the same (L, v )-sector. 

In the state (|3TiD, the degree of the pre-exponential polynomial, L, which is the total 
angular momentum, is redistributed between internal excitations, J=L—v, and the contri- 
bution due to the collective motion, v, which we call seniority for brevity. 

The energies of the states (|30|) , i.e., the eigenvalues of @ have the form 

£ k (L,v)=N + L + E h (L,v) (32) 

where Ek(L,v), the interaction energy, comes from diagonalization of V. It is therefore 
meaningful to consider the ground state as a function of both L and v, as illustrated Fig. 3. 
We call these states, 

V (L,v), £ (L,v) = min k {£ k (L,v)} (33) 

11 



with minimum energy at fixed values of L and v the "generalized yrast states" p2[ , J23 . 



The usual yrast states are those among the (|33|), which minimize the energy Sq(L,v) with 
respect to the seniority, 

£yrast{L) = min v {£ (L,v)}. (34) 

At given L, one has exactly L allowed values of the seniority, they are 

v = 0,1,2,...,L-2,L. (35) 

The value v = L — 1 is excluded because the symmetric degree J = 1 polynomial in variables 
Zi — Z is evidently reduced to zero, 

N 

polyl(zi - Z) = J2(Zi - Z) = 0. 

The states with definite seniority of type (|3CiD , can be obtained by applying the seniority 
projector [Cf. Eq.flgg)] 

V v = — d(j)e i<t>{ - v - A+A > (36) 
2ir Jo 

to the states of type flU]). 

The quantum number seniority helps to establish very useful "corresponding rules" be- 
tween the states in different (L, v )-sectors and to relate their spectra. Indeed, by virtue of 
and we can write 

V^ k (L + l,v + l) = VA + ^ k (L, v) = 
= A + VV k (L, v) = E k {L, v)A + V k (L, v). (37) 

Eq.(p7[) means the following: let {^ k (L,v)} be the p{L) normalized eigenstates for the 
Hamiltonian (0) in the sector L. Then exactly p(L) eigenstates in the sector L + 1 (having 
nonzero v !) can be obtained simply by applying the collective ladder operator to the states 
in L sector, 

/ v \ 1/2 

tt fe (L + l,v + l) = (-— - J A + V k (L,v), (38) 
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while their interaction energies will be the same, 

E k (L+l,v + l) = E k (L,v). (39) 

Similar relationships appear in the fermionic problem of the fractional quantum Hall effect 
PI- 

As the total angular momentum grows from L — 1 to L, the only new structures in the 
wave functions appear in the seniority zero sector, where we have new states of total 
number equal to 

g (L)=p(L)-p(L-l), (40) 

They must be obtained by diagonalization of V in the sector v = 0. 

Using (|40|) recursively from v to v + 1, we calculate the total numbers of states, g v (L), 
in the (L, w)-sectors 

g v (L) = p(L — v) — p(L — v — 1) for v < L — 2 

and q v= l{L) = 1. (41) 

These important relations will be very usefull in the analysis given in the following sections. 



V. ALGEBRAIC DECOMPOSITION OF THE HAMILTONIAN 



Regular methods to obtain exact ground state without solving the whole spectrum are not 



available. We use the approach [24] which we loosely nicknamed "algebraic decomposition" 



pq| . Suppose that the Hamiltonian can be written as a sum, 

V = V + V s , (42) 

such that 

(a) the first term, Vo, is simple, and one can find out its ground state |0) with eigenvalue 
£mim possibly degenerate. If the second term, Vs, has the two properties: 
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(b) Vs annihilates the state |0), so V5 10) = 0, 

(c) Vs is non-negative definite, Vs > 0, (it does not have negative eigenvalues), 

then the state |0) will still be the ground state for the full Hamiltonian Vo + Vs, with the 
same eigenvalue £ m i n - Indeed, (b) implies that |0) is an eigenstate for the sum Vo + Vs 
with its eigenvalue intact. As one knows from linear algebra, if an Hermitean operator Vq 
is perturbed by a non-negative definite Hermitean operator Vs, the eigenvalues can only 
increase (see, e.g. f29j). This means that the states other than |0) can gain energy [Cf. (b)]. 



As |0) is already the ground state for Vq, it will be the same for Vq + Vs |30 



A trivial but helpful example is the one-dimensional linear oscillator with the Hamiltonian 
V = 1/2 + a + a. The state |0), which obeys a|0) = 0, is "the ground state" for V = 1/2 
because any other state \n) has the same eigenvalue 1/2. The non-negative definite operator 
Vs = a + a, Vs > can be considered as "perturbation". It removes this degeneracy and 



leaves |0) with lowest eigenvalue of the total Hamiltonian V pO 



The same arguments apply to the case with additional conserved quantum number v, 
such as [v, Vq] = [v, Vs] = 0. In the above scheme, the single state |0) is replaced by the set 
of states \L,v), each having the property (a) and (b) in their w-sectors. This is illustrated 
in Fig. 4, where the spectrum of Vq is taken degenerate in each v-sector. 



VI. ALGEBRAIC TRIAD 

The starting strategy is to look for states annihilated by a part of V which can be then 
shown to be non- negative definite operator. In our case, the "algebraic triad" Vo, V5 and 
|0) can be established by inspecting action of terms V& in (B7|) on selected states of partition 
basis fll4|) . We observe first, that using the algebraic identity 

N 1 N N N 

N at at = - J2(at - a+)(a; - a d ) +^ a tYl a v 

i i,j i j 

one can cast the first term, Vq, of the operator expansion ( p7|) in the following form, 

N 

Vo = -j[(N - l)s — (L — A + A)( Sl - S2 )]. (43) 
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We recall now that the operator A + A is diagonal in the seniority basis (p0[). In this basis, the 
operator V is therefore reduced to a combination of quantum numbers, fl43|), with A + A = v. 
Next, one can see that the simplest basis state, \L), [Cf. Eq. (|HD1 with partition [1, 1, ...1] [] 

\L) = P s a+a+...a+\0), (44) 

is annihilated by the remainder of the Hamiltonian (j27D , 

{V-V )\L) = 0. (45) 

In order to prove (f|^), we use Eqs. (|7|,|23D to evaluate the commutator 

[a? 2 ,a+a+] = -2(a+a 12 + 1). (46) 

From this equation, it follows that 

a k 12 afatat...at\0) = 
for any k > 2, by virtue of (§). Therefore, we have 

B k \L) = (47) 
for any k > 2. From the same commutator fl4*6|), we have 

[ai2 a i2i a t a t]\°) = [ a i2 a i2, at a t]\°) = - a i2 2 |°)> 
and thus (B 2 — 5 X )|L) = 0. This relation together with (17) results in (15). 



Equations ( ^3|) and (|45f) hint that the state |L) and the operator Vo may constitute 
important ingredients of the algebraic decomposition, Cf. Eq.(42|) described in the previous 



section. However, the state \L) does not have definite seniority and it can not be eigenstate 
of the Hamiltonian. Nevertheless, as we will see, the \L) turns to be a generating function 



for states with definite seniority. Indeed, substitution Zi = Z{ + Z [see (30Q] transforms \L) 
to a sum 



x An example of such states is given by the first line of (|l5|). 
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P s ~z l ~z 2 ...~z L + ... + Z L - 2 P s ~zi~z 2 + Z 



|0> 



E W> (48) 

i;=0 



of exactly L states of the form fl30|), each being the eigenvector of A + A and therefore of Vo, 
with the eigenvalue fl45|). We notice that each t> -sector is represented by a single term in 
(0), identified with V V \L). 

The spectrum of Vq is very simple. A schematic example is shown in Fig.4, left hand side. 
It consists of L equidistant (except v^L—1), g v {L)-io\d degenerate levels with energies given 
by (|43l). The degeneracies are given by Eq.([|l]). Each v-level contains one and only one 
state \L, v) = V V \L) from the sum (|48|) . Therefore, the set of states \L, v) obey the criterion 



(a) with the operator V . The property (b) with Vs=V—V holds by virtue of (|45|). In 
particular, (a) together with (b) mean that \L, v) are the eigenvectors of V=Vq+Vs with 



eigenvalues fl43|) . This holds for any interaction V(r). 

The "algebraic decomposition" with the triad V , Vs and \L, v) = V V \L) would be com- 
plete if we succeeded to prove non-negative definiteness Vs>0 of the remainder of the Hamil- 
tonian V — Vq, criterion (c). So far, we did not specify form of the interaction V(r) in (P?D. 
We will study now general case and specify the class of potentials which have Vs > 0. 



VII. NON-NEGATIVE DEFINITENESS OF "PERTURBATION" 



We have to check signs of all the eigenvalues of Vs in the partition space ([14]). To avoid 
solving the whole spectrum in the space of symmetrized states, we use the following trick. 



By definition, the nonzero eigenvalues of Vs in the space (14) coincide with the nonzero 



eigenvalues of PsVsPs in the full space of monomials m in (14). This latter space 
has dimensionality much higher than p(L) and it includes wave functions of all possible 
symmetries, including boson sector (fully symmetric), fermion sector (fully antisymmetric) 
etc. In this extended space, the analysis of signs of eigenvalues is however crucially simplified, 
while the contributions from the symmetric sector can be accurately separated. Using the 
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symmetry of V$ under permutations of particles, we can write 

N(N — 1) 



P S V S P S = PsY, V s,U p s 



PsVs^Ps 



(49) 



where Vs,ij is the contribution from pair of particles i, j to Vs [Cf. ( pT7p6] , f42"| )1. In order to 
see that 



PsVsPs > 



(50) 



for a given interaction potential V(r), it is sufficient to show that Vs t i 2 >0, because applica- 
tion of any projector Ps from both sides in fl49|) can add new zero eigenvalues, but can not 
add negative eigenvalues. This follows from the known "inertia theorem" of linear algebra 
| 29fl . Now, we study the eigenvalues of Vs,n- Let 7r 12 be the operator of permutation of 
variables 1 and 2. The triad 

T = {n 12 , af 2 a 12 , V s ,u} 

forms a set of mutually commuting operators. Indeed, Vs,iz is expressed in terms of B>\ 2 = 
af 2 [see (pTDj. For fixed pair of particle indices, the operators a\ 2 and a\ 2 commute like 
Bose operators, [012,012] = 1. Consequently, any operator B>\ 2 can be expressed in terms of 
af 2 ai 2 , using the standard boson calculus formula 



a 12 a 12 = ai2 a 12( a l2 a 12 - l)---(a] t 2 a 12 ~ k + 1). 



(51) 



The triad T is diagonalized simultaneously in the basis of monomials Tfi — Z^ Z<2 Z^ . . . Z jsj j 
(0) with the only substitutions z\— > -^{z\—z 2 \ z 2 ^^(zi+z 2 ), 



V2 



(zi - z 2 ) 



h r 



V2 



Zi + z 2 ) 



4 3 -^|o). 



(52) 



In this basis, the eigenvalues of the triad T depend only on l\ through the subfactor 



ll of the eigenvector (|52|) . These eigenvalues are given by 
T = {(-1) 1 \ h, X h }, 



(53) 
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respectively, where A; 1 is the eigenvalue of Vs,i2 which can be readily calculated. Using 



Eqs.([27|), (|26|), (|51"|) and the summation formula 



we obtain the expression for the eigenvalue of Vs^n in the form of the integral 

oo 

K = I rdrV(r)f h (r) (54) 

with / the functions defined as 



o 



1 + 4(2- — 



(55) 



2H X \ 4 V 4 

The eigenvalue (— l)' 1 of 7r 12 helps now to separate out the states with wrong symmetry: the 
eigenvectors with li odd are antisymmetric in Zi,Z2, and the projector P$ in fl4~9|) eliminates 
their contributions. All even values of l\ (<L) can contribute to the bosonic sector, and 
the corresponding A's must be checked. Therefore, the set of inequalities for the control 
eigenvalues 

A 2n > (56) 



for all values of integer n obeying 2n<L forms the sufficient condition for (|50|), i.e., for the 
"perturbation" Vs to be non-negative definite, criterion (c) Under these conditions, the 
triad V , Vs = V — V and {V V \L)} obeys the "algebraic decomposition" with properties 
(a), (b) and (c), with V V \L) being the ground state in its sector L, v with the energy S min 
equal to eigenvalue of V . 

The conditions ( |56"|) can be used to describe the universality class of the interactions 
V{r). They are analyzed explicitly in the following sections. 



2 Note that f (r) = and f 2 (r) = 0, therefore A = A 2 = 0. 
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VIII. GENERALIZED YRAST STATES AND THEIR SPECTRA 



With the criterion (|56|) and Eqs.(2B) and (f43|) at hand, we can formulate very general 



result: For any bona fide two-body potential V(r) which satisfies the integral condition 



CO 



t 2n / f T 

1 + n 1 - 



> 



(2n)\ \ 2 

for any n < Ll /2, (57) 



the eigenstates of the Hamiltonian H with minimal energies at given pair of v and 



L(< min{Z/, N}) have universal form, as follows from (^8 

/ a \N-L+v n 

1 N 

2-]y5>- (58) 

Their energies are given by (|43D and they can be expressed through the simple moments of 
V(r), 

S min (L, V ) = L + N+ ^'^ so + '^^Niv - L) (59) 



s 



roo i>oo I \ +2\ 

/ dte-'ViV^t), S! = / dte- 1 - + - V(V2t) 
Jo Jo \2 A J 



which are equal to expectation values of interaction between two bosons both in the ground 
s o = (V'oV'ol^l'^oV'o) and the first excited state of oscillator si = > respectively. 

Here, (^^{V^ip) denote the usual two-body matrix element and ipo = |0) and ipi = z\0). It 
is worth to note that this relation does not mean that the eigenvalues ( (5"9"D can be obtained as 
an expectation value over a simple state, say, with the two lowest oscillator levels occupied; in 
fact, the eigenstates ([58]) are much more complicated when written in the second quantization 
representation. 

At fixed L, we have exactly L such equidistant generalized yrast states, marked by v = 
0, 1, 2, L (v ^ L — 1), Cf. Fig. 3., right hand part. Each such state is the "ground state" 
in the sector L, v (of course, there are other states in each sector with higher energies). 
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Example of the spectrum of a real system for N = L = 6 and the (^-interaction is shown 
in Fig. 5. 

The usual yrast states are those of fl5!f ) that minimize S m i n (L, v) with respect to v, Cf. 
Eq.fl3ip. As £ m i n (L, v) in fl59|) depend linearly on v, it is immediately seen that there is only 



two options f] : 

(I) v = 0, if D = s - si > 0, (60) 

(II) t» = L i/ D = s -si<0. (61) 
where the "spectral discriminant", D, is given by 

D = jT°° 'dfe-*fi-^y(V2^. (62) 

For the 6ona /ide potentials obeying (|5^), the first option is usually realized, as is assumed 
in Figs. 3 and 4, see also Figs. 5 and 7 for particular interaction potentials. For example, 
integrating by parts in (|62"D, the spectral discriminant D is seen positive for any decreasing 
potential 

^1 < (63) 
dr 

which allows the representation The condition (|63D can be replaced by 

dV(r) 



dr 



< 0, (64) 



if strict inequality in ( |63"D holds for at least single value of r. 

The inequality (Vo — Vi)(v — L) < corresponding to (I) means that internal rotational 
excitations with higher J = L — v are energetically favorable, once the interaction energy 
between two bosons in the state z\0) is smaller than that in the state |0). 



3 The rare case of exact equality in ( |6l| ) corresponds to complete degeneracy of the generalized 
yrast states with different values of v, having the same N and L, see, e.g. Ref. [p4|). 
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Physically, the yrast wave functions (pq) with v = correspond to condensation to a 
vortex, rotating around the "center-of-mass" , as discussed, for example, in ||. 

In contrast, the maximum seniority states with v = L, which correspond to purely col- 
lective rotation Q with no internal excitations, were shown [[| to be energetically favorable 
in the case of attractive <5-forces. It is curios that there exists a broad class of predominantly 
attractive interactions, for which the yrast states have the same form as for the attractive 
5-function interaction |2f| . 

The results desribed in this section are very general |22] . The condition ([57]) defines the 
class of the potentials for which the results (|58|j59|) are valid. In the next section, we show 
that (|57D holds for many potentials of physical interest and consider few explicit examples. 
Discussion of general properties of ( pT[) will be given in section 10. 



IX. THE RESULTS FOR PARTICULAR INTERACTIONS 



We consider now few applications of the results obtained in general form in the previous 
section to particular cases of various potentials of interaction between the bosonic atoms. 
Gaussian forces with variable range and S -interaction. 
We start with considering the case of repulsive Gaussian interaction, 

U 



V(r) 



7TR 2 



U >0, 



(65) 



with Uq a nonnegative strength and R the radius which can be varied from zero to infinity. 
We obtain from Eq. ( [57] ) the control eigenvalues 



A 



Un 



2n 



tt(2 + R 2 ) 



R 2 



2/i 



- 1 + 4n- 



1 + R 2 



(66) 



2 + R 2 J (2 + R 2 ) 2 

From this equation, it can be easily seen that the non-negative definiteness condition (|5?D 
^2n > is fulfilled for any n irrespectively to the value of R, see Fig. 6. This can be easily 
proved by induction in n. Indeed, we have from (|66|) taht Aq = A2 = and 



A 



2(n+l) 



\2n 



4U (1 + R 2 
tt(2 + R 2 ) 



R 2 
2 + R 2 



2n 



> 0, 
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therefore \2 n > for any n > 2. The results given by Eqs.( [58| , |5"9"|) are therefore valid for 
any L < N, and the spectrum of the generalized yrast states for the case of the Gaussian 
potential is given by 



n(2 + R 2 ) 



N(N — l)/2 



S min {L,v) = L + N + 

;N(L — v) 



(67) 



(2 + i? 2 ) 2 

One can see that the usual yrast states correspond to v = 0. 

Taking the zero range limit, R — > 0, in Eqs. ( |65|,|66 , |67 ), we pass to the case of the 
5-function repulsive interaction 

S(r) 



V = U 5(r) = U 



2nr 



(68) 



In this case, we obtain instead of (|6 



Uo 
2vr 



^1,0 + 



(69) 



It is seen that while Ai=— 1/2, for any l\=2n even the condition \<i n > (57) holds, Cf. Fig. 
6. The energies of the generalized yrast states ( |5§D are now given by 

(70) 



S mm (L, v)=L + N+ ^N(2N — L + v — 2). 

8n 



Yrast states have v = 0, and £ m i n (L, 0) agrees with that obtained numerically [f|, see also 

0, a, m- 

It is expedient to look at the full spectrum of the system. At low iV and L, the problem 
can be diagonalized analytically. Below, we present the results for the total angular momen- 
tum L = 6 in the system of six particles interacting via repulsive 5-function interaction. The 
total number of states ( dimensionality of the partition basis ) is p(6) = 11. The interaction 
energies Ek(L,v) of the states are given by 



E (6, 6) = 15—, Eq(6, 4) = 12—, E (6, 3) = — — , 



^o(6,2) =9 



Uo 
2vr' 



^(6,2) 



21Uo 
2 2tt' 
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2 2vr' v ' ' 8 2vr : 



£ (6,0) = ^6, ^(6,0) = ^ 



69 ; 

+ IT + C 



£ 2 (6,0) = | 



69 > 



with 



£ 21 

32 + 4£ ; 



«M)-£(* + f) 



(71) 



£ = (324 + 12iv / 32199) 1 



/3 



Here, the index marks degree of excitation of the state within the same L, w-sector: so 
Eq(L,v) = E min (L,v) correspond to the generalized yrast states (|59| , |70D . There are no 
excited states in the sectors v = 6, 4 and 3, the sectors v = 4 and 5 have one excited state 
each. The sector v = has three excited states, their energies are found by solving a cubic 
equation. The spectrum is shown in Fig. 5. 
Two- and three dimensional Coulomb forces 

Of special physical interest are the long-range forces of Coulomb type. For the conven- 
tional Coulomb interaction 



V 



Ua 



U >0 



(72) 



we obtain in ( |57j ) 



A 



2ra — <-^0 



C/oV2tt 



'bn 1 r(2n + l/2) > 



(73) 



16 2 20F(2n)! J ' 

It is easily seen, using induction in n, that the inequalities A 2n > Q57| ) hold (Cf. Fig. 6), 
and the wave functions of the generalized yrast states are therefore given by (pSQ. Their 
energies are given by 



S min (L,v) = L + N + 



The yrast states correspond to v — 0. 



U V2tt 
64 



N[16(N - 1) - 5{L - v)}. 



(74) 
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Similar formulas can be obtained for the two-dimensional Coulomb log-interaction 



V = U, 



Jog Q 



U >0, 



which corresponds to repulsion at small distances. We have 



A - Uo 

*2n — —T 



3n + 2^(1) - 2^(2n) 

n 



> 0, 



(75) 



(76) 



where ip is digamma function p7|, the inequality in \2 n > for ( fflf) can again be easily 
proven by induction, Cf. Fig. 6, so condition fl57D holds. The energies of the generalized 
yrast states are given by 



L + N 



U N 



(/^(2)-7)(JV-l)-|(i;-L) 



(77) 



with 7 = — = 0.57721... the Euler constant. The yrast states correspond to v = 0. 

The results very similar to the above can be obtained for the screened Coulomb (Yukawa) 
forces V oc exp(r/r )/r. 

The results for the generalized yrast states in the cases considered in this section are 
illustrated in Fig. 7 for N = L = 6. 



X. THE SUFFICIENCY CONDITION AND THE UNIVERSALITY CLASS OF 
BONA FIDE REPULSIVE POTENTIALS 



The simple sufficiency condition (57) can be checked straightforwardly for any interaction 
potential V(r) of interest, as was done in the previous section. It is interesting to understand 
why many potentials do meet this condition. In this section, we show that the applicability 
condition ( |57D imposes in fact only weak restrictions on the class of forces V(r) which can 
be regarded as predominantly repulsive. We will discuss the condition ( 1571 ) in more details 
and define some subclasses of the bona fide potentials of interest. 

Short-range interactions 
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The factor functions f2n(r) in expression J °° rdrf2 n (r)U(r) of Eq. ([57]) are plotted in 
Fig. 8. At small r ~ 0, the factor-functions approach positive constant values with zero 
derivatives, 

f n \ =n-l, f^l =0, (78) 

r=0 df r=0 

while the first node of the functions /2n( r ) occurs at 



r = V2V / 12-2v / 30~ 1.43 (79) 

which is of order of the oscillator length, — , in our units (h — m — uj — 1). From Eqs.(|T8"D 
and ([75]), it follows therefore that for any short range (as compared to the characteristic 
length of the trap) interaction, the condition ( |5"7| ) reduces to the single inequality 

J d 2 rV(r) > 0, (80) 

Thus, the results ( |58| , |5"9"|) hold for short-range interactions, which are repulsive on average. 
I is seen also that inequality (j30|) gives v=0 for the usual yrast states, according to (|60|) . 

It is interesting that the sort-range potentials do not have to be purely repulsive to match 
the condition (|57|), (|80|) . Instead, the condition ([57]) implies that V(r) has sufficiently strong 
repulsive component. For example, consider the two-parameter family of the potentials in 
the form 

V{r) = -^-{l-{r/Rf)e-f^ (81) 

These potentials can be regarded as attractive in the usual sense, having the potential well 
at short distance. Examples are shown on Fig. 9. In fact, if the well is deep enough, it can 
even support bound state. Indeed, consider the Schroedinger equation of relative motion 
(s-wave) of the pair of particle with unit mass each, interaction via fl5I|), 

r dr \ dr J 

for the wave function of the ground state, 0o- We obtain rigorous variational upper bound 
on the ground state energy, e , using the trial wave function O = exp(—^r 2 /R 2 ): 
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V6Vo \720V* + 1 - 8(6V ) 5/2 /3 - 46V - 2(6V ) 3/2 
e < l - 



R 2 (2AV - l) 2 

It is seen that we have negative energy of the ground state, Eq < 0, if 



216 



V > 

1 2 



(3763 + 18V502) 1/3 + 241 , , + 13 

(3763 + 18y/502y/ 3 



~ 8.99 (82) 

in such cases the potential well indeed supports the bound state(s). On the other hand, the 
potential ( |8~TD satisfies the basic condition (0) if 



R < — V5 + 5V41 ~ 0.608, (83) 

with all its control \ 2n nonnegative, irrespectively of the magnitude of Vo, Cf. Fig. 9, 
lower panel. Therefore, the family of the potentials (|8T| ) satisfying ([52]) and ( p3|) will be 
"predominantly repulsive" in the sense of the condition (|57D, having the solutions for the 
generalized yrast states in the form (p8|) and (|59|) , while having strong attractive component 
sufficient for the bound state. In the sense of three-dimensional scattering theory, the 
scattering length can be negative. 
Long-range forces 

The condition ( |57|) holds even for many long-range interactions. This is seen from be- 
havior of function fi{r) (n = and n = 1 give / = 0) which is positive at r < ro and 



r>ri = y2V 12 + 2^30 ~ 3.10, (84) 

and /4 is negative only in the interval ro < r < r\ (f n (^) for higher n behave similarly). By 
direct calculation, it is easy to see that the integrals of the functions rf 2n [Cf. Fig. 8, lower 
panel] over r vanish, 

oo 

/ rdrf 2n (r) =S 1 + S 2 + S 3 = 0, 
o 

ro n oo 

Si = J rdrf 2n (r), S 2 = J rdrf 2n (r), S 3 = J rdrf 2n (r), (85) 



r rx 
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so the net areas coming form the regions where / > and / < coincide, 

Si + S3 — IS2I) 

as is shown in Fig. 10, upper panel, for the case f±. From this geometry, which is also 
illustrated in Fig. 10, lower panel, it is clear that fl57|) holds, if V(r) decreases monotonically 
and fast enough, as is the case for the long range Gaussian and Coulomb forces. Example 
of the potentials of such type are shown in Fig. 10, lower panel. 

Integrating by parts, one sees that fl57D holds for any interaction described by monoton- 
ically decreasing and concave function V(r) 0, 

not necessarily repulsive everywhere. This means, in particular, that the solutions fl57|) apply 
to the exponential and screened Coulomb potentials, 

V e (r) = \V \exp(-r/R), V sc {r) = \V \ ^^- . (87) 

One can therefore summarize that the condition holds for any physically meaningful 
repulsive interaction. 

XI. CONCLUSION 

To conclude, we considered the problem of weakly interacting Bose atoms in symmetric 
harmonic trap. Our main purpose was to study the yrast states of the system. In order 
to treat the problem of the ground state of the system at a given angular momentum, we 
developed an expansion of the interaction in powers of ladder operators. This universal 
operator expansion is exact and it provides a very convenient way to study the eigenvalue 
problem in the coordinate representation. 



*It is assumed that the potential V can be represented in the form of operator expansion (E7|). 
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Taking into account the additional conserved quantity, associated with the collective 
contribution to the total angular momentum, we considered a more general problem, namely, 
the ground state of the system as a function of two conserved quantum numbers: the total 
angular momentum and the angular momentum of internal excitations. We called this series 
of states "generalized yrast states" . 

A method of "algebraic decomposition" for the interaction has been developed to derive 
analytically the states with minimum energy at a given angular momentum and seniority 
(generalized yrast states). The results apply, in particular, to the usual yrast states. The 
wave functions of the generalized yrast state have the form of "condensed vortex states". 
Their energies are expressed through their quantum numbers and simple integrals of the 
interaction potential. 

We studied the condition on the interaction potential which allows the use of our solution. 
Analysis shows that there exists a broad universality class of the repulsive interactions for 
which these results are valid. We described this universality class by simple integral condition 
on the interaction potential. 

The results of the work allow further generalizations and developments. The three- 
dimensional case can be done using the same method. It is also interesting to study region 
of higher angular momenta L > N, where the structure of the basis polynomials will be 
changed f26|, while the numerical studies indicate signs of phase transition ||. The method 
of "algebraic decomposition" developed here is not restricted to this particular problem and 
can be applied to fermions and even to the particles with parastatistics. 

The work was supported by FAPESP. 
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Figure Captions 



Fig. 1. Definition of the yrast states. The spectrum of a system (circles) is shown versus 
the angular momentum, L. The lowest energy states at each value of L, connected by the 
curve, compose the "yrast line" of the system. 

Fig. 2. Fragment of the spectrum of the system (schematic plot) within approximation 
of the weak coupling limit. The left hand side is the spectrum of the system without 
interactions. The levels are equidistant with the spacing %u, each level is p-fold degenerate. 
Central part is the spectrum of interacting system. Each L-level splits onto p levels. The 
sequence of the states with minimum energy at given L is the yrast-line (right hand part). 

Fig. 3. Illustration to the definition of the generalized yrast states. The left-hand side 
is the part of the spectrum of the interacting system with definite angular momentum L 
(corresponding to a bunch of levels in the central part of Fig. 1.). These states can be 
grouped onto "bands" with definite values of the seniority, v , as is shown on the right hand 
side. The sequence of lowest energy states in their w-sector composes the generalized yrast 
states (they are marked with asterisks). One of those, with absolute minimum of energy, is 
the usual ground state. 

Fig. 4. Illustration of "algebraic decomposition" . The spectrum of Vq (left) is sequence 
of degenerate levels, labeled by the conserved quantum number v. The non- negative definite 
perturbation V$>0 splits each level, pushing the states up and leaving the lowest energy in 
each v-sector intact. The resulting spectrum of Vo + Vs is shown on the right. 

Fig. 5. The spectrum of the system with repulsive 5-interaction calculated analytically 
for the case of six particles in the sector with total angular momentum L = 6. The total 
number of states is p(6) = 11. The energy of the levels e(L,v) = ^[S(L,v) — L — N], 
Eq. flT0|) , are plotted against the internal angular momentum, J = L — v. The generalized 
yrast states drop on the straight line described by Eq.([7D|), they are connected by dashed 
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line. 



Fig. 6. "Redge trajectories" for A/ x (0). Upper panel: the values of a = ^A; are plotted 
versus / (circles connected by curves) for two different interaction potentials. The curves 
correspond to / considered as continuous variable. Solid curve - repulsive Gaussian potential 
with R = 1 and dashed curve - repulsive ^-function interaction. The double circles denote 
the control eigenvalues with even I = 2n, Eqs.( |57| . pE| ). 



Lower panel: the same but for the Coulomb interaction (solid curve) and Zog-Coulomb 
interaction (dashed curve), Cf. Eqs.( [73]J7B| ). 

Fig. 7. The specta of the generalized yrast states e = j^S min (L,v) are plotted against 
J = L — v (N = L = 6) for the Gaussian interaction with R = 1 (small circles connected by 
solid line), 5-interaction (small circles, dashed line), Coulomb interaction (big circles, solid 
line) and the log-Coulomb interaction (big circles, dashed line). 

Fig. 8. Upper panel: Factor-functions /2„(r) for n = 2,3,4, 5 are plotted against r : / 4 
(curve 1), f 6 (curve 2), / 8 (curve 3) and f w (curve 4). 
Lower panel: the same as above but for the r/2 n (r). 

Fig. 9. Upper panel: the potentials [Eq.(jSl~l)1 supporting the bound state with Vo = 10 
obeing Eq . (|81~1) are plotted against r for three values of R: R = 0.3 - solid curve, R = 0.6 - 
dashed curve, and R = 1.0 - dashed-dotted curve. 

Lower panel: the "Redge trajectories" for (0). The eigenvalues A/ in combinations 
a = j^A; are plotted as function of I (circles) for the potentials shown on the upper panel. 
The curves connecting the symbols correspond to eigenvalues as continuous function of I: 
R = 0.3 - solid curve, R = 0.6 - dashed curve, and R = 1.0 - dashed-dotted curve. The 
double circles denote the control eigenvalues with even / = 2n. 



Fig. 10. Upper panel: under curve r/ 4 (r) (illustration to Eq. 

Lower panel: the function r/ 4 (r) (thin curve) is plotted against r together with potentials 
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matching condition (|57|). The latter are shown by thick curves: 1 - Gaussian potential 
V(r) = ^exp(-r 2 /R 2 ) with R = U = 1.5, 2 - Coulomb potential V(r) = ^ with 
U = 1/5 and 3 - log-Coulomb potential V(r) = Uolog(l/r) with U = 1/3. 
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Fig. 2. Hussein & Vorov "Condensed... " 
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Fig. 3. Hussein & Vorov "Condensed..." 
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figurei Hussein & Vorov Fig. 4. "Condensed... ' 
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Fig. 5. Hussein & Vorov "Condensed..." 
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Fig. 7. Hussein & Vorov "Condensed..." 




Fig. 10. Hussein & Vorov "Condensed ..." 




